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1.  Introduction 


Within  the  armor  development  community,  there  is  a  need  to  develop  a  compact, 
disposable,  high-fidelity  blast  gauge  to  measure  overpressures  associated  with 
harsh  events,  such  as  deflagration  of  energetics  or  detonation  of  high  explosive 
where  accelerated  projectiles  and  debris  may  occur.  Many  times,  overpressures 
generated  by  such  events  can  be  a  nuisance  to  neighboring  systems  and 
environments.  Knowledge  of  the  pressure  history  can  allow  engineers  to  design 
systems  to  mitigate  these  effects.  There  are  numerous  gauges  and  techniques  that 
have  been  developed  to  measure  blast  pressures;  however,  these  systems  are  often 
expensive,  too  bulky,  too  complicated,  or  too  inaccurate  to  use  in  desired 
implementations.  For  instance,  current  piezoelectric-based  blast  pressure  pencil 
probes  from  PCB  Piezotronics1-3  provide  a  solution,  but  the  gauges  cost  near 
$2000,  are  constructed  of  352  g  of  metal  (plus  a  BNC  cable  connection),  are  greater 
than  300  mm  in  length,  and  only  work  within  specified  ranges  of  overpressures 
necessitating  some  amount  of  prior  knowledge  of  the  expected  overpressure  to  be 
measured. 

This  report  describes  the  construction  and  functionality  of  a  Photonic  Doppler 
Velocimetry-  (PDV-)4,5  based  blast  pressure  gauge,  which  applies  enhancements 
within  the  field  of  velocity  measuring  technology  to  improve  blast  measurement 
technology.  The  gauge  is  relatively  cheap  costing  near  $50,  is  constructed  of  mostly 
soft  plastics  and  thin-glass-based  fiber  optic  connectors,  can  be  scaled  in  size 
depending  upon  the  desired  measurement,  and  is  only  limited  in  measurement  range 
by  the  yield  strength  of  the  disk  material  on  which  the  gauge  is  based. 

2.  Results  and  Discussion 


The  PDV-based  pressure  gauge  works  on  the  principle  that  a  blast  wave 
propagating  through  a  medium  such  as  air,  incident  normal  to  a  nondefonnable 
cylindrical  disk,  will  accelerate  the  disk  following  Newton’s  second  law  of  motion. 
For  strong  overpressures,  Peng  et  al.6  identified  a  computational  method  capable  of 
determining  the  blast  pressure  that  accounts  for  the  pressure  increase  generated  by 
a  shock  at  the  disk’s  rear  surface.  This  situation  is  pictorially  described  in  Fig.  1.  If 
the  disk  mass  and  shape  are  known  and  unchanging,  and  complications  associated 
with  release  waves  generated  near  edges  of  a  viable  disk  of  finite  dimensions  are 
minimized,  the  acceleration  history  of  the  disk  is  described  by 

dud  _  Apd  .. 

dt  pdhd  '  '' 
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Fig.  1  Schematic  depicting  blast  pressure  interactions  with  a  free-flying  plate 

Where  m  is  the  disk  velocity,  t  is  time,  Apd  is  the  difference  in  pressure  at  the  front 
and  rear  surface  of  the  disk,  pd  is  the  plate  density,  and  hd  is  the  plate  thickness.  By 
comparing  simulations  of  the  disk  response  to  the  experimental  measurements  of 
the  disk’s  acceleration,  one  could  back  out  the  blast  pressure  history  applied  to  the 
puck.  In  reality,  the  analysis  would  necessitate  employment  of  a  complete 
3 -dimensional  geometry  describing  the  local  environment,  or  many  simulations  of 
finite  time  steps,  as  the  blast  pressure  history  may  not  follow  a  simplistic  linear  or 
exponential  form  as  simplified  in  Peng  and  colleagues’  computational  validation.6 

Figure  2  shows  a  schematic  (A)  and  a  photograph  (B)  of  the  device  designed  to 
make  the  aforementioned  measurement.  The  device  primarily  consists  of  a  metal 
cylinder  suspended  at  the  opening  of  a  polymer  tube  using  3  skinny  posts.  A  PDV 
probe  is  mounted  at  the  opposite  end,  which  is  used  to  measure  the  acceleration  of 
the  disk  as  it  is  accelerated  from  the  blast  pressure.  This  devise  was  constructed  to 
accomplish  the  following  constraints: 

•  The  disk  should  not  yield  elastically  during  the  event. 

•  The  device  should  be  capable  of  making  an  accurate 
position/velocity/acceleration  measurement  of  the  disk  throughout  the 
event. 

•  The  accelerations  of  the  disk  from  nonblast  events  should  be  minimized. 
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Fig.  2  Schematic  (A)  and  photograph  (B)  of  a  blast  pressure  gauge  that  uses  a  heterodyne 
velocimetry  measuring  technique 

The  first  constraint  is  dependent  upon  the  intensity  of  the  blast  being  measured.  For 
relatively  low-pressure  fields,  such  as  that  generated  by  release  of  compressed  air 
from  a  standard  shop  air  compressor,  the  disk  could  be  constructed  using  most 
materials  such  as  plastics,  woods,  metals,  and  so  forth.  For  higher  intensity  blasts, 
such  as  that  generated  by  release  of  energy  from  a  high-explosive  detonation  or 
deflagration,  materials  such  as  metals  or  ceramics  may  be  needed.  A  relatively 
robust  initial  disk  design  applicable  for  blast  pressures  similar  to  that  generated  by 
detonation  of  high  explosive  of  up  to  10  lbs,  TNT  equivalent,  at  an  offset  distance 
of  1  m  uses  a  cylindrical  disk  machined  from  mild  steel  of  radius  6.4  mm,  thickness 
2.05  mm,  and  weight  2.0  g. 

To  address  the  second  constraint,  a  single  PDV  channel  was  used  to  measure  the 
acceleration  history  of  the  disk  along  the  devices  symmetry  axis.  This  system 
allowed  for  a  high- temporal  profile  (~1  ns)  and  high- accuracy  (~  1  %-2%) 
measurement  of  the  disk.  Because  PDV  is  limited  in  its  identification  of  the  source 
versus  target  movement,  it  is  necessary  to  ensure  that  the  PDV  probe  is  held 
stationary  during  any  measurement  event,  or  that  its  motion  is  negligible  compared 
to  the  translation  of  the  nondeformable  disk  when  acted  on  by  the  blast  front.  In 
general,  probe  motion  can  be  reduced  by  increasing  the  inertia  of  the  probe  itself. 
This  can  be  accomplished  via  sturdy  coupling  of  the  PDV  probe  to  the  blast 
measuring  device  and  attaching  to  a  large  mass  via  dispersive  media  such  as  a  foam, 
neoprene  rubber,  or  Isodamp.  Figure  3  shows  one  example  of  how  such  a  device 
was  coupled  to  a  large-mass  metallic  frame  when  assessing  the  blast  loading  of  a 
Multi-Energy  Flash  Computed  Tomography  system. 
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Fig.  3  Depiction  of  a  robustly  mounted  blast  pressure  gauge,  with  a  shock  isolating  Neoprene 
layer,  used  to  measure  shock  loading  on  an  aluminum  support  structure 

Using  a  single  PDV  probe  enables  the  possibility  of  error  if  the  disk  is  rotated  or 
translated  off  of  the  symmetry  axis  by  a  nonradially  unifonn  pressure  field.  If  this 
were  the  case,  however,  the  simplistic  relation  to  the  pressure  field  described  in 
Eq.  1  would  be  invalid  and  a  more  complicated  analysis  would  be  necessary.  One 
could  be  alerted  to  this  situation  by  collecting  multiple  co-timed  PDV 
measurements  from  a  few  radial  locations  on  the  disk,  or  using  a  full  field 
diagnostic  (i.e.,  optical  photography)  to  record  the  disk  movement. 

To  address  the  third  constraint,  efforts  were  made  to  ensure  that  the  disk  was 
initially  mounted  in  an  axis-symmetric  fashion,  the  disk  was  weakly  coupled  to 
initial  external  mounting,  and  influences  that  could  perturb  a  planar  or  spherically 
symmetric  blast  front  were  minimized  throughout  the  measurement.  Of  these 
efforts,  ensuring  that  the  disk  was  weakly  coupled  mounted  in  a  planar,  axis- 
symmetric  fashion  was  found  to  be  the  most  important.  By  design,  the  gauge  was 
constructed  using  3  threaded  rods  to  “press-fit-hold”  the  disk  in  space.  The  threaded 
rods  allowed  for  both  tension  adjustment  of  the  “press-fit”  and  precise  control  to 
center  the  disk.  To  assess  the  best  method  of  coupling  the  disk  inside  the  holder, 
empirical  observations  of  the  disk  release  were  made  using  high-speed  photography 
for  many  methods  of  disk/threaded  rod  coupling  when  subject  to  an  air  blast 
generated  using  an  air  compressor  and  nozzle.  This  data  is  summarized  in  the 
following  Table  and  shown  in  Fig.  4.  The  coupling  methods  included  1)  threaded 
rods  in  as-is  form  (not  shown),  2)  threaded  rods  prepared  with  conical  points  to 
minimize  contact  with  the  disk  (not  shown),  3)  using  an  off-take  of  Whisker  Biscuit 
bow  and  arrow  rest  in  which  fiber  brushes  were  used  to  initially  hold  the  disk  in 
place,  4)  solid  grooved  ring,  5)  adding  flexible  rubber  caps  to  the  ends  of  the  rods, 
6)  adding  springs  to  the  ends  of  the  threaded  rods,  and  7)  adding  spring-coupled 
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ball  bearings  to  the  ends  of  the  rods  to  reduce  the  friction  while  sustaining  a  suitable 
coupling  force  prior  to  the  release  from  mounting.  Of  these,  the  spring-coupled  ball 
bearing  holder  appeared  to  have  the  cleanest  release,  although  the  air  loading  was 
not  symmetric,  causing  the  disk  to  translate  and  rotate  in  an  nonaxisymmetric 
motion.  For  the  particular  configurations  tested,  6-mm  nylon  or  steel  ball  bearings 
were  epoxied  to  the  ends  of  3-mm-outside-diameter  springs  constructed  of  0.5-mm 
steel  wire,  which  were  capable  of  exerting  a  maximum  of  1.16  kg  of  force. 

Table  Description  of  the  disk  holding  methods  and  their  qualitative  results 


Method 


Qualitative  result 
(10=good,  l=poor) 


Threaded  rods  3 

Conical  pointed  rods  4 

Brushes  . 

(“Whisker  Biscuit”) 

Solid  w/grooved  ring  5 

Rubber  caps  7 

Springs  1 

Springs  w/ball  bearings  9 


Description  of  result 

Rigid  rods  were  inflexible  resulting  in  compressive 
loading  along  the  rod  length  as  the  disk  attempted  to 
release.  In  some  cases  this  precluded  the  disk  from 
releasing. 

Only  slightly  better  than  non-modified  threaded  rods 
resulting  from  smaller  surface  contact  with  the  disk. 
Disk  had  increased  drag  during  release  as  it  “worked”  its 
way  through  the  bristles. 

Clean  release,  but  complicated  by  possible  jetting  of  gas 
through  the  opening  of  the  ring. 

Clean  release. 

Flexibility  of  springs  resulted  in  inability  to  firmly  and 
flatly  hold  the  disk  prior  to  shot. 

Very  clean  release.  Difficult  to  load  and  keep  symmetric. 
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Brushes  ("Whisker  Biscuit®") 


Solid  w/Grooved  Ring 


Springs  w/  Ball  Bearings 


Fig.  4  Photographs  of  numerous  techniques  used  to  loosely  hold  a  disk  within  the  blast 
pressure  gauge  and  the  images  depicting  how  the  disk  was  ejected  when  subject  to  an  air  blast 

To  reduce  device  influences  on  the  blast  front  and  its  interactions  with  the 
nondeformable  disk  throughout  measurement,  efforts  were  made  to  minimize 
materials  interacting  with  the  blast  front  in  the  vicinity  of  the  disk  (with  the 
exception  of  the  solid-grooved  ring  holder).  In  particular,  the  3  threaded  rods  that 
hold  the  disk  were  symmetrically  spaced  and  minimized  in  diameter.  We  can  also 
delay  influences  of  reflected  waves  by  maximizing  the  boundary  standoff  of  the 
diagnostic  body  (tube)  from  the  disk.  This  is  pictorially  described  in  Fig.  5  where 
Sv  is  the  wave  front  velocity. 
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Pressure 

Front 


Fig.  5  Schematic  depicting  the  shock  propagation  geometry  within  a  blast  pressure  gauge 


The  distance  a  reflected  scattered  wave  must  traverse  (xrefiected)  is  represented  by 
the  tightly  dashed  black  line  and 

xrefiected  =  / Sin  d)  ~  2  r  ^ '/ Sin  o)  ■ 

Where  d  and  Q  are  defined  in  Fig. 5,  T,  is  the  tube  radius  and  Pr  is  the  disk  radius. 
The  distance  the  disk  travels  (xd)  in  this  time  is  proportional  to  the  integral  of  the 
integral  of  the  force  applied  by  the  blast  front  (F)  divided  by  the  mass  of  the  disk 
0 m ): 

xd°c  f*oV(u)du,  (3) 


V  (it)  oc  f  t  —  du  .  (4) 

Although  we  cannot  solve  these  equations  without  additional  information  about  the 
specific  impulse  and  materials  involved,  we  do  observe  that  T,~Pr  is  in  the 
numerator  of  Eq.  2  and  the  disk  mass  is  in  the  denominator  of  Eq.  4.  Therefore,  if 
we  either  increase  the  ratio  of  tube  radius  to  the  disk  radius,  or  decrease  the  disk 
mass,  the  distance  over  which  a  pure  measurement  can  be  made  will  be  extended. 

To  assess  the  validity  of  the  PDV-based  pressure  gauge  design,  an  experiment  was 
perfonned  in  which  1 14  g  of  Primasheet  1000  was  detonated  adjacent  to  a  gauge 
that  used  the  spring/ball-bearing  coupling  method.  The  gauge  was  positioned  0.8  m 
from  the  center  of  the  blast  as  shown  in  Fig.  6.  The  high-explosive  charge  was 
constructed  of  multiple  sheets  stacked  to  make  an  almost  rectangular  geometry,  and 
was  detonated  from  the  top  using  a  Teledyne  RP-80  explosive  bridge  wire 
detonator.7  From  bottom  to  top  the  high-explosive  stack  included  5  pieces  of 
Primasheet  1000  C-5  50  mm  x  50  mm,  3  pieces  of  Primasheet  1000  C-l  50  mm  x 
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50  mm,  and  one  piece  of  Primasheet  1000  C-l  25  mm  x  50  mm.  The  total  stack 
thickness  of  the  50  mm  x  50  mm  sheets  measured  to  29.7  mm,  and  the  total 
thickness  including  all  high  explosive  measured  30.5  mm. 


Fig.  6  Photograph  of  the  experimental  setup  used  to  validate  a  blast  pressure  gauge 

The  PDV  signal  derived  from  Doppler-shifted  light  reflected  off  of  the  rear  surface 
of  the  disk  was  acquired  using  a  Third  Millennium  Engineering8  FI  77a  mod  block 
coupled  with  a  16-GHz  Keysight9  digitizer  bandwidth  limited  to  1.25  GHz.  The 
signal  was  analyzed  using  MATLAB-based  Sandia  InfraRed  HEtrodyne  aNalysis 
(SIRHEN)  software  package,10  which  uses  a  Fourier-based  frequency  analysis 
method.  This  combination  of  system  and  analysis  provided  capability  to  measure 
up  to  approximately  1  km/s  with  a  temporal  resolution  of  2*1  O'6  s. 

Figure  7  shows  a  multipanel  figure  including  (A)  the  PDV  spectrogram,  (B)  the 
pertinent  region  of  disk  acceleration,  and  (C)  calculation  of  the  difference  in 
pressure  fields  applied  to  the  disk.  Figure  8  overlays  the  acceleration  of  the  disk  as 
measured  with  that  simulated  using  Velodyne.11  Figure  9  displays  extracted 
pressure  fields  from  tracers  imbedded  within  the  simulations  5  mm  in  front  of  the 
disk,  at  the  disk  surface,  and  10  mm  behind  the  disk  surface.  For  the  remainder  of 
the  analysis,  we  focus  on  the  movement  of  the  disk  in  the  region  outlined  by  the  red 
dashed  box  overlaid  on  Fig.  7A.  During  this  portion  of  the  measurement,  the  disk 
translated  approximately  1  mm  from  its  initial  position. 
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Fig.  7  Multipanel  image  showing  (A)  the  spectrogram  resultant  from  the  blast  pressure  gauge 
from  the  experiment  shown  in  Fig.  6,  (B)  the  extracted  acceleration  calculated  by  taking  the 
derivative  of  data  extracted  in  panel  A,  and  (C)  the  computed  AP  to  which  the  disk  within  the 
blast  gauge  is  subjected 
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Fig.  8  Overlay  of  the  experimentally  measured  disk  acceleration  and  the  computed 
(Velodyne)  disk  acceleration  (red)  and  probe  acceleration  (blue) 


Fig.  9  The  pressure  fields  to  which  the  disk  within  the  blast  pressure  gauge  was  subject  as 
computed  by  Velodyne 
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Experimentally,  the  disk  begins  accelerating  approximately  490  ps  after  detonation 
of  the  high  explosive,  which  computes  to  a  linear  approximation  of  the  shock 
velocity  traveling  through  air  at  1.63  km/s.  This  is  significantly  greater  than  the 
elastic  compression  velocity  in  air  at  standard  temperature  and  pressure 
(0.343  km/s).  The  disk  reaches  a  peak  velocity  of  6.36  m/s  approximately  200  ps 
after  its  initial  acceleration.  During  this  period,  it  was  found  to  undergo 
accelerations  up  to  57000  m/s2  (5800  g’s),  which  corresponds  with  a  difference  in 
pressure  near  9*10A5  Pa  (130  PSI).  It  is  worth  noting  that  accelerations  of  this 
magnitude  are  well  resolved  by  the  capabilities  of  the  PDV  system. 

Comparisons  of  the  simulations  with  the  experimental  data  demonstrate  similar 
features  within  the  disk  acceleration  traces:  a  smooth  acceleration  uptake  of 
matching  slope,  general  overall  shape  to  include  slight  general  deceleration 
(presumably  from  drag  and  turbulence  associated  with  the  flow  of  high-explosive 
gases  in  air  and  reflection  of  expanding  gasses  off  of  the  gauge  wall  surfaces),  and 
similar  temporal  duration  of  features.  The  simulations  reveal  that  the  gauge 
geometry  allows  the  pressure  front  to  interact  with  the  rear  mount  that  holds  the 
PDV  probe  approximately  120  ps  after  the  disk  initially  moves  (blue  trace  in 
Fig.  8).  Because  PDV  was  fielded  in  the  conventional  manner,4  one  is  not  able  to 
discern  experimentally  if  velocity  changes  in  the  spectrogram  are  a  result  of  the 
reflector  accelerations  or  the  probe  accelerations.  From  this  analysis,  one  can  infer 
that  the  large  deceleration  recorded  in  the  experimental  data  near  610  ps  is  a  result 
of  the  probe  being  influenced  by  the  pressure  front. 

The  absolute  magnitude  of  the  simulated  disk  acceleration  peaks  at  83%  of  that 
measured  with  the  PDV  blast  gauge.  The  PDV  technique  is  known  to  be  capable  of 
measuring  a  reflector’s  velocity  within  l%-2%  of  its  absolute  value,  indicating  that 
simulation  shortages  likely  account  for  the  remaining  15%  error.  Because  the  disk 
acceleration  has  a  linear  dependence  on  the  difference  in  pressure  fields  as 
described  in  Eq.  1,  one  can  infer  that  the  pressure  fields  are  underestimated  by 
approximately  15%.  This  discrepancy  may  be  due  to  the  approach  used  in  Velodyne 
to  model  fluids  and  fluid  structure  interaction.  For  example,  using  a  peak 
acceleration  of  35000  m/s2,  one  would  expect  a  corresponding  pressure  of  550  kPa. 
However,  the  same  simulation  predicts  a  maximum  peak  pressure  of  350  kPa.  This 
suggests  an  issue  with  how  the  pressure  influences  the  motion  of  a  small,  relatively 
thin  solid  body. 

Applying  the  aforementioned  factors  to  the  pressure  field  simulations  displayed  in 
Fig.  9  suggest  that  the  pressure  field  generated  by  the  experimental  geometry  peaks 
in  the  range  of  8-9.2*  105  Pa  upon  interaction  with  the  disk  surface.  This  value  is 
only  sustained  momentarily  as  the  pressure  field  is  increased  by  a  factor  of  2  from 
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reflection  at  the  surface,  but  dissipates  quickly  back  to  a  lower  value  near 
6-6.9*  105  Pa  as  rarefactions  propagate.  These  values  compare  well  with  the 
experimentally  assessed  9*  1 05  Pa,  as  the  shock  generated  behind  the  disk  surface 
in  this  geometry  is  expected  to  be  small. 

3.  Conclusions 


A  PDV-based  blast  gauge  was  designed  that  provides  a  cost-effective,  high- 
resolution  solution  toward  measuring  overpressures  generated  during  blast  events 
where  accelerated  ejecta  and  debris  may  be  an  issue.  The  gauge  benefits  from  recent 
enhancements  in  measuring  free  surface  velocities  using  heterodyne  PDV-based 
techniques.  The  gauge  design  allows  for  low  impact  on  experiments,  applicability 
to  a  wide  range  of  pressures,  scalability  in  size,  and  is  disposable. 

Within  the  gauge  design,  optimization  was  perfonned  to  find  the  most  suitable 
coupling  method  to  hold  the  disk  in  place  prior  to  interaction  with  an  overpressure 
wave.  It  was  found  that  a  spring/ball-bearing  coupling  method  produced  the 
cleanest  release  of  the  disk  with  the  least  amount  of  induced  rotation. 

The  gauge  using  the  spring/ball-bearing  coupling  method  was  fielded  in  a  high- 
explosive  detonation  experiment,  demonstrating  the  ability  of  the  gauge  to  measure 
blast  overpressure.  In  comparisons  to  simulations,  the  response  mechanics  of  the 
gauge  compared  well;  however  the  simulations  under  predicted  the  accelerations 
(and  pressure)  by  approximately  15%.  These  values  may  be  improved  if  codes  such 
as  ALEGRA12  or  ALE3D13  were  used  to  simulate  the  event. 
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